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Abstract. We investigate perturbed Hamiltonian systems with two degrees of 
freedom by renormalization group method, which derives a reduced equation called 
renormalization group equation (RGE) by handling secular terms. We found that 
RGE is not always a Hamiltonian system. The necessary and sufficient condition 
that RGE becomes a Hamiltonian system up to the second leading order of a small 
parameter is that the original system is separable by Cartesian coordinates. Moreover, 
RGE is intcgrable when it is a Hamiltonian system. These statements are partial 
generalizations of our previous paper. 



Short title: Canonical structure of renormalization group equations 
February 5, 2008 
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1. Introduction 



To understand the nature of non-linearity in dynamical systems, it is indispensable to 
investigate temporal evolutions of the system. Particularly, in Hamiltonian systems, 
slow relaxation of 1/f spectra is observed in systems with two |T], Q and many || |[] 
degrees of freedom, and understanding its cause is an important subject of Hamiltonian 
dynamical systems. The slow relaxation is connected with self-similar hierarchical 
structure of phase space, which is observed in a 2-dimensional mapping ||. The 
structure is also supposed in systems with many degrees of freedom with the aid of 
Kolmogorov- Arnold- Moser theorem || and Poincare-Birkhoff theorem ]7| . By assuming 
self-similar hierarchical structure and transition probabilities between levels of the 
hierarchy, some abstract models || || |K| are proposed to understand the slow relaxation. 
The abstract models catch universal properties of Hamiltonian systems, but they are 
not available to investigate details of individual systems owing to their assumptions. 
For instance, Poincare recurrences P(t) (integrated probability to return into a given 
region after a time larger than t) decays with t as P(t) ~ l/t p , and the exponent p 
takes different values in different 2-dimensional mappings 0. One approach to grasp 
such individual properties is to numerically calculate equations of motion, and another 
approach is to analytically construct approximate solutions. In this paper we focus on 
the latter approach in nearly integrable systems with perturbation techniques. 

Naive perturbation often yields secular terms due to resonances and they invalidate 
perturbative expansion. This problem of secular terms is dealt with singular 
perturbation techniques fTIf , e.g. averaging methods, multiple scale methods and 
matched asymptotic expansions. However, these techniques are not easy to use 
because we must select a suitable assumption about the structure of perturbation 
series. Recently, perturbative renormalization group method |12[ [13| is proposed as 
a tool for global asymptotic analysis of the solution of differential equations. It unifies 
the methods listed above, and can treat many systems irrespective of their features. 
Renormalization group method omits fast motion and constructs reduced equations 
of motion from secular terms. The reduced equation is called renormalization group 
equation (RGE), which is geometrically regarded as an envelope equation. 

Although RGE gives approximate but global solutions to original systems, we 
have one question: whether does RGE keep properties of original systems? We miss 
characteristic behaviors of orbits in phase space if RGE becomes a dissipative system 
which often has attractors and orbits exponentially converge to them. In contrast, 
Hamiltonian systems never have attractors and orbits eternally wander. 

In the previous paper [14|, we treated a Hamiltonian system with a homogeneous 
cubic or quartic potential function added as perturbation, and showed that if the RGE 
becomes a Hamiltonian system then the original Hamiltonian system is integrable. By 
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investigating properties of RGE, we can justify integrability of the original Hamiltonian 
system. The present paper aim to extend the result for general perturbation function 
V(qi, #2) • Moreover, we mention integrability and an integral of RGE. 

The plan of this paper is as follows. In section |^ we review renormalization group 
method by using a simple example. The main theorem is shown in section ^, and it is 
proved in sections £|, || and [| Section |7| is devoted to summary and discussions. 



2. Renormalization Group Method 

Here we review renormalization group method by using a simple system. Let us consider 
the system 

H(q,p) = H {q,p)+eV(q), (1) 

where 

Ho(q,p) = l(p 2 + q 2 ), V{q) = ^q\ (2) 
and e is a small parameter. The equation of motion is 

and the exact solution is 



q = B cos(VTTe t) + C sm(vT+e t), (4) 

where B and Co are constants of integration and determined by initial condition at the 
initial time t = to- 

We perturbatively solve equation (g) by expanding q as a series of positive powers 

of e: 

q = q + eq x + t 2 q 2 H . (5) 

This naive expansion gives 



q(t; t , B , C ) = B cost + C sint + -(t — t ) (Co cost — B sint) 



€ 



2 



(t — t ) (C cost — B sint) + (t — t ) 2 (-Bo c °st + C sint) 

o L 
+ 0(6 3 ). 



(6) 



We ignored other homogeneous parts of qj (j > 1), which are kernel of the linear operator 
L = d 2 /dt 2 + 1, since they can be included in go- This naive expansion breaks when 
e(t — t ) > 1 because amplitude of eq\ exceeds one of q Q . 
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To remove the secular terms, we regard B and C$ as functions of the initial time to. 
The temporal evolutions of B(t ) and C(t ) are determined by the following equation 

rani: 



dq 
df 



= 0, for H. (7) 

t =t 



From the equation (@), we obtain the RGE 
dB(t) (e 



dt 
dC{t) 



C(t) + 0(e d ), 



e e 



2 N 



(8) 



- ---\B(t) + 0(e 



3^ 



dt V 2 

where -B(t ) = -Bo an d C(^o) — Co- The renormalized solution q RG is 
q RG = q(t;t,B(t),C(t)) 
= B cos 




and this solution is the same as the exact solution ([|) up to 0(e 2 ). 

Let us check whether the reduced equation (|8|) reflects properties of the original 
system (ffl), which are (i) symplectic properties, (ii) integrability and (iii) the form of an 
integral. In equation (H), B and C are canonical conjugate variables and the equation 
is governed by the Hamiltonian 

(10) 




and canonical equations 

dB(t) _ dH RG 

~ (11) 
dC{t) dH RG V ; 

dt ~ dB~' 

Equation (|^) is therefore integrable and an integral of the system is B 2 + C 2 which 
corresponds to H Q . As shown in section |5| these statements hold in any Hamiltonian 
systems with one degree of freedom whose integrable part is a harmonic oscillator. What 
happens in systems with two degrees of freedom where chaos may appear ? This is the 
main question of this article and we show a part of the answer in the next section. 



3. Main theorem 



In this section we present the main theorem. 
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Theorem 1 Let Hamiltonian be represented as follows: 

H(qi,q2,pi,p 2 ) = H (q 1 ,q 2 ,pi,P2) + eV(q u q 2 ), (12) 

where 

H (q l ,q 2 ,p 1 ,p 2 ) = ^{pl+pl + ql + ql), (13) 

e is a small parameter and the potential V(qi, q 2 ) is an analytic function which has only 
even order terms of qi and q 2 , for instance, q\q\. We perturbatively solve the equation 
of motion by setting the 0-th order solution as 

q j0 (t) = Bj cos t + Cj sin t. (j = 1, 2) (14) 

Then the following two conditions (A) and (B) are equivalent. 

(A) The renormalization group equation (RGE) of the system ({Hj) is a Hamiltonian 
system up to the second order of e. 

(B) The original system ( fUj j is separable by Cartesian coordinates, that is, rotation of 
the coordinates q\ and q 2 . 

Moreover, when (A) (or (B)) is satisfied, the RGE is integrable and (Bj+Cf+B%+C%)/2 
is an integral which corresponds to H . 



Remark 1: RGE for the system ( |T2| ) is always a Hamiltonian system up to 0{e l ). Hamil- 
tonian of the RGE is the same as one obtained by canonical perturbation theory [RJ 
and the Hamiltonian is time-average of V{q\o, q 2 o) as shown in subsection |4.2| . 

Remark 2: RGE is an equation for variables which are constants of integration 
introduced in the 0-th order solution, and hence the form of RGE changes when we 
use different form of the 0-th order solution, for instance, 

q j0 (t) = ^2AjCos(t + 9j). (15) 

However, the theorem holds for new variables if they are canonical conjugate pair, since 
constructing RGE and changing constants of integration are commutable operations. 

Remark 3: We excluded odd order terms of q± and q 2 from the perturbation V(q\,q 2 ). 
We must construct RGE up to 0(e 4 ) if V(qi,q 2 ) has odd order terms (c.f. [0) and 
need longer calculation. But the extension is straightforward. 

We prove this theorem 1 in sections |], |5] and |(J Sections [| and [5] are devoted to 
prove (A)=^(B) and (B)=^>(A), respectively. Section ^| is for integrability and the form 
of an integral. 
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4. Proof of (A) (B) 

We first derive RGE up to 0(e 2 ) for system (03), and we investigate conditions that RGE 
is a Hamiltonian system. Then we compare the conditions with separability conditions 
of system flT2|). 



4-1. Renormalization Group Equation 
Equation of motion in system (|I2|) is 
qi + qi = -eV 1 (q 1 ,q 2 ), 

(16) 

h. + ?2 = -£^2(91,92), 

where the subscripts of V represent partial derivatives with the variables, namely, 
Vi = dV/dqi and Vi 2 = d 2 V/dqidq 2 - We expand (ft and g 2 as power series of e, 

gi = gio + egn + e 2 gi 2 + 

^2 = <?20 + e<?2i + e 2 ?22 H • 

Equations of motion in each order of e are 

Q30 + Qjo = 0, (18) 

Qji + Qji = - Vjtfo), (19) 
2 

qj2 + qj2 = - J2 v jk(%)qku (20) 

k=l 

where go = (?iO)92o) and j = 1,2. Following an assumption described in theorem 1, we 
set the 0-th order solution as 

q j0 = Bj cos t + Cj sin t. (j = 1, 2) (21) 

Since Vj(q ) is a periodic function of time with period 2n, we expand it by Fourier series: 

oo 

- Vj (g ) = } + { u t ] cos rat + v j n) sin • ( 22 ) 

n=l 

Then solution to 0(e x ) is 

(0) . (* — *o) / (1) . , (1) ,\ . 1 / (n) , . (n) . ,x /oq n 

g^i = ; H — (it] sint- t>J- cos tj + } j [Uj cos nt + sin nt) , (23) 

2 n=2 1 — n 

where to is the initial time. 

Let us proceed to solution to 0(e 2 ). We also expand Vjfc(go) by Fourier series, 

oo 

- V jk (q ) = u$ + cosnt + sinnt ). (24) 

n=l 



Terms in the right-hand-side of equation ( p0|) which make secular terms in qj 2 are 

2 

— E Vjk(qo)Qki — Gjt cos t + Hjt sin £ + K 3 - cos t + Jj sin £, 

k=l 



where, by using relations between Fourier components shown in Appendix A 





= 2£ 




k=l 




2 


Hj 


= 2E 




k=l 




d 


K 3 






dB 3 




d 







/ Q U (0) Q2 U (0) ^ Q U (0) Q2 U (0) 



( du^ <9V°) _ du<® d 2 u 



(0) 



i u 



(0)\2 



(4 0) ) 2 



[11 



n=2 L 
oo 2 



n- 



{ut ] ) 2 + (v[ n) ) 2 + (4 n) ) 2 + (v[ n) ) 2 



4 U} ) 2 + (4 Uj ) 2 + E [(^S n) ) 2 + (^ n) ) 2 + (^ n) ) 2 + (v^) 2 

, v „=2 1 n 

Solution to 0(e 2 ) is written with these symbols as 
q j2 = (t - t ) 

+ (t - to) 2 









+ — ij sint 


[(1 


-y) cost + 








V 4 





4 



cos £ H — r- sin £ 



4 



+ (non-secular terms). 



From the secular terms of equations ( P3|) and (|30|) , we obtain RGE as 

,2 



3 2 1 \ 4 2 

J 2 J V 4 2 



4 



<9$ 



as, ■ 4 



31 



where the potential $ is defined by 

,2 C oo -i 

A \ n=2 1 n 

and is the constant component of Fourier series of — V{%) 

1 r 27T 



(25) 

(26) 
(27) 
(28) 
(29) 



(30) 



(31) 



(™)\2l 



,(32) 



u 



(0) 



2tt 



V(q 10 , q 2 o) d£. 



(33) 



4-2. Conditions that RGE is a Hamiltonian system 



We show the conditions that the RGE ( |3"1"D becomes a Hamiltonian system. The 
conditions are 



Ai = A 2 = A 3 = A 4 = A 5 = A 6 = 0, 



(34) 
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where (/ = 1, 2, 



Ai = ^ + 



dB 1 



6) are defined by 
A 2 = 



dB 2 



dB 2 dC 2 



A, 



dB 2 
dBi 



+ 



<9C 2 ' 
dC 2 



A s = 



A, 



dB~2 



+ 



<9B 2 



(35) 



dC 2 dBi 



dC 2 dCi ° dB 2 3d 
The term written by $ in equation always satisfies conditions (|33j). Accordingly, 
whether RGE is a Hamiltonian system is determined by Gj and Hj, which originates 
in secular terms of the 1-st order solution. From equations 
forms of A^'s become 



and 



the concrete 



Ai 



A 5 



m 

2 2 / 



q u (o) q du (o) q 
V 9B k dC k 
du<® d 



d 2 



k=l 



A, = A 4 = — 



V dB k dC k 



dC k 3B k/ 

dvF> d s 
dC k dB k/ 

( d 2 



dB 2 

' d 2 

dB 2 



+ 



d 2 



dc\ i 
d 2 ' 



u 



(0) 



■it 



(0) 



d 2 



0B 1 dB 2 \8C 2 



u 



(o) 



2 / 



d 2 



d 2 



2 fdu® d 9u(°) 9 9uW d 



2 ( <9.B a dd dB x 



+ 



dddC 2 \dB 2 
du<® d \ 



dB 2 



u 



(0) 



A' 



2 



8B l dB 2 



+ 



U 2 



dB 2 dC 2 
' d 2 
dB x dB 2 
/ 92 



+ 



u 



(0) . 



dC 2 dB 2 
d 2 \ ' 
dddC 2j 
d 2 



(36) 



u 



(0) 



dddC 2 



\dB l dC l 

d 2 d 2 



+ 



0B 2 dC 2 

d 2 



u 



(0) 



d 2 ' 
9C? 



if 



(0) 



where 



A 5 
A fi 



A' + A' 



6 • 



A5-A6 



(37) 



Note that RGE is always a Hamiltonian system up to 0{e 1 ) since terms of 0(e l ) do 
not appear in Aj's. In this order, a Hamiltonian of the RGE becomes time-average of 
V(q ): 



RG 



-U 



(0) 



1 

2^ 



2tt 



dt V(q 



(3f 



and if]^ G is the same as the one obtained by canonical perturbation theory ||15, [16 



^.5. Condition that the original system is separable 

In this subsection we clarify the condition that the original system ([12]) is separable by 
rotation of the coordinates (qi,q 2 )- The integrable part H is always separable by any 
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rotation, and hence we only have to consider the perturbation V. In the next subsection, 
we compare the condition written by with equation ([36]). 



Theorem 2 The condition that the original system ( \H\ ) becomes separable by rotation 
of the coordinates (51,92) is 

Ci(V n -V 22 )+2c 2 V 12 = 0. (c 2 + c 2 = l) (39) 

In particular, when the potential V(qi,q 2 ) has only even order terms of q\ and q 2 , the 
condition $3l\ ) is equivalent to each of the following three conditions: 

d 2 d 2 \ „ / d 2 d 2 ' 



(V 

(W 

(III) 



' ( d 2 d 2 

C\ I 77777 + 



Cl 



dBl dC\ dB\ dC 2 2i 

d 2 o 2 d 2 d 2 ' 
+ 



dB\ dCl dB 2 dC 2 2i 

' d 2 d 2 N 



2c 2 
2c 2 



dB 1 dB 2 + dC 1 dC 2> 
' d 2 d 2 " 



dB l dC l dB 2 dC 2 



+ c 2 



dB 1 dB 2 dddC 2 
d 2 d 2 1 



m (0) = 0, (40) 
w (0) = 0, (41) 



+ 



dB x dC 2 dC x dB 2 



u 



(°) = 0. 



(42) 



Corollary 3 From (I) and (II), the two conditions 



(IF) 



Cl 



Cl 



a 2 


d 2 \ 


+ 2c 2 


( d 2 


ae 2 


dB 2 2 ) 






d 2 \ 


+ 2c 2 


( d 2 


dC 2 


dC 2 ) 


\dC x dC 2 



«(°) = 0, 
u<°> = 0, 



(43) 
(44) 



are also equivalent to each of (I), (II) and (III). 



Equation (|39|) is a special case of Darboux equation |]17l , which is used in the pre- 



vious paper [14]. To prove theorem 2, we introduce the following lemma. 



Lemma 4 Let us assume that V(qi, q 2 ) has only even order terms of q\, q 2 and V(qi, q\ 
is not constant. We define the function F(q 1 ,q 2 ) as 

F(Qi,Q2) = Ci\Vu(qi,q 2 ) - V 22 (gi,g 2 )] + 2c 2 V u (qi,q 2 ), 



then 



(a) F(q 1 ,q 2 ) = for y q u y q 2 

& (b) dtF(q 10 , q 20 ) = for v 5 1; V C 1; v 5 2 , V C 2 
J 

&(c) dtF(q 10 ,q 2Q ) cos2t = for y B 1 , v C 1 , V B 2 , V C 2 
Jo 

p2tt 

<*(d) / dtF(q l0 ,q 20 )sin2t = for v B u v C u v B 2 , y C 2 
Jo 



(45) 

(46) 
(47) 

(48) 

(49) 
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Proof of Lemma 4: 

The direction (a) =>- (b),(c),(d) is obvious, and hence we prove the reverse direction, 
that is, (b),(c),(d) =>- (a). We may set C\ = C 2 = since (b),(c) and (d) are satisfied 
for any values of Bj and Cj. Then F(q w , q 20 ) is expanded as 

F(q 10 , g 20 ) = E E rr 2 *!*! ™* m+n t, (so) 



where 



am+n rp 

F, """ >(0 - ) = «(°- 0) - < 51) 



As a result, in case of (b), 

0= F(q w ,q 20 )dt 
Jo 



^ F( m '™)(0, 0) nm nn \ 27i2- m - n {m+n) C {m+n)/2 (m + n: even) 



" i^o m!n! 1 2 \ (m + n : odd) 1 } 

and F^ m '^(0, 0) = for even m + n. Consequently, F(qi, q 2 ) = because F(q 1 , q 2 ) has 
only even order terms of q% and q 2 . We can prove (c), (d) =^ (a) through similar ways. □ 

Now let us prove theorem 2. 

Proof of Theorem 2: Let us prove the first half of theorem 2. We rotate coordinates 
from (qi,q 2 ) to (qi,q 2 ) as 

qi = qi cos 8 + q 2 sin 9, 
g2 = — <?i sin 6 + q 2 cos 6*. 



The separability condition is that there exists a real value of # such that 

0, forW (54) 



d 2 V 



dqidq 2 

This condition is equivalent to 

(Vn-V 22 )sm20 + 2F 12 cos20 = O, forVV (55) 
and consequently, we have proved the first half of theorem 2. 



From relations between Fourier components shown in |Appendix A| , the conditions 
(I), (II) and (III) are equivalent to 

(I): dtF(q w ,q 20 ) = 0, (56) 

Jo 
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(II) : / ^ dtF(q l0 ,q 20 ) C os2t = 0, (57) 

Jo 

(III) : r dtF(q w ,q 20 ) cost sin t= - f " dtF(g 10 , g 20 ) sin 2t = 0, (58) 

Jo 2 jo 

respectively. Then, from lemma 4, we conclude that theorem 2 has been proved. □ 
4-4- Final step of proof of (A (B) 

We prove (A)=^(B) of theorem 1 with the aid of theorem 2 and corollary 3. The 
conditions (F) and (IF) are rewritten in matrix form 

_2L_ _aLw(o) o^l 




-dB\ dB\' U ^aB,flRn U \ ( C\ \ ( 



(59) 



The constants c\ and c 2 are not zero simultaneously, and the conditions (?) and (II' 
must be degenerate accordingly. The degenerate condition is 

d 2 



det 




d 2 



8B x dB 2 \dCl dCi dddC 2 \dB 2 dB. 



u {0) = 0, (60) 



and this separable condition is equivalent to A 3 = A 4 = 0. We also find that the 
degenerate condition between (I) and (III) is equivalent to A 6 = 0. Consequently, 

RGE is a Hamiltonian system up to 0(e 2 ) 
Ai = A 2 = A 3 = A 4 = A 5 = A 6 = 
=>- Potential function V(qi,q 2 ) is separable by Cartesian coordinates. 



5. Proof of (B) (A) 

We assume that the Hamiltonian system (0) is separable by Cartesian coordinates. 
RGE for the Hamiltonian system is also separable since constructing RGE and rotating 
coordinates are commutable operations. As a result, what we have to prove is the fol- 
lowing theorem. 

Theorem 5 RGE up to 0(e 2 ) for the system with one degree of freedom 

H(qi,Pi) = ^pl + ql)+eV( qi ) (61) 
is a Hamiltonian system. V is an analytic function of q\ . 
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Proof of Theorem 5: RGE for system (|BID is obtained from equation fl3~l|) by omitting 
terms with subscript 2: 



Cl -~m +e T' 



(62) 



where 



2 ( oo i 

$ = - «*W - y |(-S 0) ) 2 + E if^[(«i n> ) a + (^ n) ) 2 ]} > ( 63 ) 

/ du (0) Q2 U (0) du^d 2 U^ \ 

Gl - 2 { dB 1 dB l dc l ~ dd dB\ ) ' (64) 

1 ~ 2 {dB 1 dC\ dC x dB 1 dC 1 J- (65) 
The condition that the RGE fl62|) becomes a Hamiltonian system is 

dB x dC 1 . . 

^ s m + *k = a ' (66 » 

and the form of Ai is 



e 2 9Ci 



2 



9 / d 2 <9 2 \ (0) _ (cP_ d_ 

dB x dCt [dBf + dCl) U ' ~ ~dC 1 ~dB~ l \dB~f + dC\ 

~0B 1 ~~dC 1 ~ ~ dd dB 1 ' 



(67) 



We prove that the Ai is always zero by showing both and u[i are functions of 
B\ + C\. From the definition of u^°\ 



u<°> = -- 
2tt 



1 /" 27r 

— / dtViBi cost + Ci sin t) 

!7T JO 

- / " dtVUfWl + C\ cos(t - 5)) 
lir Jo v 



2tt Jo 

1 00 yWfrn . /-27T (68) 

= -E i ^#( 5 i + ^ 2 ) nx2 " n2 ^ 

n=0 

where tan 5 = C\jB\ and V^ fc ) = d k V/dq k . Consequently, is a function of 5 2 + C\. 
We can also prove that uf^ is a function of B\ + C 2 through the same way. □ 
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6. Integral of renormalization group equation 

From the consideration of the system with one degrees of freedom, RGE with two degrees 
of freedom is always integrable when it is a Hamiltonian system since it is separable. 
Moreover, H of system (|T^) is an integral of RGE. 



Theorem 6 For the RGE system with one degrees of freedom (pzj), the quantity Bf + C( 



is an integral of the equation. For the RGE system with two degrees of freedom (\3j\), if 
RGE is a Hamiltonian system then B\ + C\ + B\ + C\ is an integral of the equation. 

Proof of Theorem 6: First we consider systems with one degree of freedom. We 
prove that the time derivative of {B\ + C 2 )/2 becomes zero: 

^V-^n^^'^ 1 < 69) 

We prove equation (|69|) by showing that $ is a function of Bf+Cf and that BiG+C\H = 
0. 

We first consider the terms written by <3>. The concrete form of $ is 

,2 f oo i 

$ = -«*<» - 6 - |4 0) + E if^K«?°) 2 + (^ n) ) 2 ]} • ( 7 °) 

The first and the second terms u^°\ uf*\ are constant components of — V(qo) and 
— Vi(q Q ), respectively. They are therefore functions of Bf+Cf accordingly (of. equations 
(p*8|)). From the definition of u^, 
1 r 2 ^ 

— I &W\(Bi cos t + Ci sin t) cos nt 

77 JO 

1 dtV^J Bf + Cf cos(t - 5)) cosnt 

7T JO 



(n) 



~ dt E ^Tp-iBl + Clf' 2 cos k {t - 5) cosnt 



1 £ Y^M {B 2 + cl f\C^ cosnS - S« sinn5), (71) 
17 k=o K - 



where 



r (k) _ 



r2n /"27T 

/ dtcos k t cosnt, S^' = / dt cos fc t sin rat. (72) 
jo jo 



We also have 



- E ^^(5? + C?)*^((7« sinn5 + cosn6) } 



71 t^o ]A 



(73) 
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and 

71 fc=OZ=0 K ' f ' 

Consequently, $ is a function of -Bj* + Cj 2 since C^C$ + S^S® is independent of £>i 
and Ci. 

Next we prove that B\G + Ciif = 0. From (64]), (|65 ) and the fact that m^ -* is a 
function of Bf + Cf, 

G = -2Ci[(?i (0) )f , F = 25 1 [(u (0) )f, (75) 

where (u^)' is derivative of u^(Bf + Cf), and hence 5 X G + Ci# = 0. 

Finally we consider systems with two degrees of freedom. When RGE is a Hamilto- 
nian system, both B\ + C\ and B\ + C\ are integrals of RGE in separated coordinate. 
Accordingly, Bf + Cf + B% + Cf is an integral in any Cartesian coordinates since it is 
invariant by rotation of the coordinates (51,52) which induces rotation of two pairs of 
variables: (B 1 ,B 2 ) and (Ci,C 2 ). □ 

Now we have completed proof of theorem 1. □ 



7. Summary and discussions 

We considered renormalization group equation (RGE) for perturbed Hamiltonian system 
with two degrees of freedom whose integrable part is two harmonic oscillators with the 
same angular frequency, and perturbation is a function of positions having only even 
order terms. RGE is not always Hamiltonian system and we presented the necessary 
and sufficient condition that RGE becomes Hamiltonian system as the theorem 1. The 
theorem states that RGE in Cartesian coordinates becomes a Hamiltonian system up to 
0(e 2 ) if and only if the original system is separable by rotation of the coordinates. When 
RGE is a Hamiltonian system, it has an integral which corresponds to the integral part 
of original system and RGE is integrable accordingly. The theorem also assert that RGE 
is a Hamiltonian system when the original system has one degree of freedom. RGE is 
always Hamiltonian system up to 0(e l ) and its Hamiltonian is the same as one obtained 
by canonical perturbation theory. 

The theorem excluded odd order terms from the perturbation potential function. 
We must construct RGE up to 0(e 4 ) to extend the theorem to any analytic perturbation 
functions since odd order terms give information on the separability at 0(e 4 ). But it is 
straightforward to extend our result to include odd order terms in potential function. 
We expect that the theorem holds for any perturbation from the result of homogeneous 



cubic perturbation case WM, which satisfies the theorem. 
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RGE is not a Hamiltonian system when the original system has two degrees of 
freedom and have no integral other than Hamiltonian of the system. On the other hand, 
RGE with one degree of freedom is always Hamiltonian system. We therefore expect 
that we can understand properties of chaotic orbits by investigating how symplectic 
properties break in the process of constructing RGE. Solutions to non-Hamiltonian 
RGE will also give some aspects of chaos. 

Future works of this topic is as follows. We suppose that it is straightforward 
to extend the main theorem of this paper to systems with many degrees of freedom. 
Another extension is generalization of integrable part of original systems. We are 
particularly interested in the case that integrable part is non-linear and angular 
frequencies change as values of actions. 
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Appendix A. Some relations between Fourier components 

We show some useful relations between Fourier components of V(qo), —Vj(qo) and 
—Vjk(qo)- From the definition, is represented by as 



7T JO 

d 
dB[ 
du<® 
dB 1 



dtVi(gi ,g 2 o) cost 
1 r 2 ^ 

- / dtV(gio, g 2 o) 

7T JO 



Similar relations hold for other components: 



- 2— 



= 2- 



and for ra > 2, 
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Some Fourier components of Vjk(qi , q 2 o) are rewritten by it™ as follows: 
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From the definition of Fourier components, the following equation is always satisfied for 
any Fourier components: 



d 2 



-X 



d 2 



-x. 



dB 1 dC 2 dC 1 dB 2 
where X is an arbitrary Fourier component. 



(A5) 
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Appendix B. Hamiltonian of RGE in quartic perturbation 

We show the explicit form of the Hamiltonian of RGE in quartic perturbation when 
RGE is Hamiltonian system. First we describe how we obtain the Hamiltonian from the 
RGE. Second the concrete form of the Hamiltonian is shown. 



Appendix B.l. Method of constructing Hamiltonian of RGE 
When RGE is a Hamiltonian system, it is written as 

h{B\, Ci, B 2 , C 2 ; e) : 



dt 
dCi 

~df 

dB 2 
dt 

dC 2 
dt 



g 1 {B l ,C 1 ,B 2 ,C 2 ;e) 
fii.Bx, d, B 2 , C 2 ] e) 
g 2 {B x ,C x ,B 2 ,C 2 -e) 



dCt ' 

dH RG 

obT' 

g H RG 
dH RG 



(B6) 
(B7) 
(B8) 
(B9) 

(BIO) 
(Bll) 

where is an arbitrary function of B\,B 2 and C 2 . To determine M\, we substitute 
equation (|B10| ) to equation (|B7|) and we obtain 



dB 2 

By integrating equation (|B6| ) with respect to C\, we have 
H RG = H RG ^ + M^(B 1 ,B 2 ,C 2 ), 



dB 1 



dB 1 



9i 



(B12) 



Ci=0 



since the left-hand-side of equation ( |B12| ) does not depend on C\ when the RGE is a 
Hamiltonian system. By integrating equation (|B12|) with respect to Bi, we get as 



where 



mV>(B 1 ,B 2 ,C 2 ) = # RG(2) + M ( - 2 \B 2 ,C 2 ), 
< g H RG(l) 



H RG(2) 



dBt 



9i 



dB±. 
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Next, to determine M^ 2 \B 2 ,C 2 ), we substitute equations (|B10| ) and (|B13| ) to 
equation (|B8|), then 



dC 9 
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(H RG ^ + H RG ^) + f 2 



(B15) 



S 1= Ci=0 
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By integrating equation (pi5| ) with respect to C 2 , we get M^> as 
M^(B 2 ,C 2 ) = H RG ^ + M^\B 2 ), 

where 

' d 



H RG(3) 



da 



. (jE RG(l) +jff RG(2) ) + /2 



dC,. 



(B16) 



(B17) 



_B 1= Ci=0 



Finally, to determine we substitute equations ( P10|) , ( |B13|) and (P16|) to 

equation (|B"8|), then 

dM^> ( d 



d£ 2 V aB 2 ; y 

By integrating equation (P18| ) with respect to B 2 , we get as 
M( 3 )(5 2 )=if RG ( 4 ), 

where 

^RG(4) = J ^__^_ (jff RG(l) + F RG(2) + jjRG^ _ \ 

Consequently, the Hamiltonian H RG is 

H RG = H RG ^ + H RG ( 2 ) + if RG ( 3 ) + if RG ( 4 ) 



S 1 =C 1 =C 2 =0 



(B18) 
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B 1= Ci=C 2 =0 



(B21) 



Appendix B.2. RGE Hamiltonian for quartic perturbation 
We present explicit form of Hamiltonian of RGE for the system 

H = \(Pi + v\ + Qi + q\) + e(axqi + a 2 q\q 2 + a-,q\q 2 2 + a^\q\ + a 5 q 2 )- 

We know that RGE for this system has Hamiltonian if if the coefficients a±, 
satisfy the following condition [IIJJ: 



(B22) 



9a 2 + 4^3 — 24oi0 3 — 9o 2 o 4 = 0, 
9«4 + 4«3 — 24«3«5 — 9a 2 «4 = 0, 
(a 2 + 0:4)0:3 — 6(040:4 + a 2 os) = ( 
Then the H RG is 



where 



H RG = Ua l {Bl + Clf + a,{Bl + a 
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+[a 2 (Bl + Cl) + o 4 (5 2 2 + C 2 2 )]{B l B 2 + C X C 2 )} 
+\a,{Bl + Cl){B 2 2 + C 2 2 ) + \oc^{B l B 2 + C X C 2 ) 



. «5 
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# 2 RG = "| [«?(*? + Clf + al{Bl + Clf] 

~ §~ 2 \^Bl + Cl)\B\ + C\) + a\{B\ + Cl){Bl + C 2 ) 2 ] 

- ^(5? + C D( B l + <3)(*i 2 + B l + ^ + Cf) 

- ^[a^Bl + C 2 ) 2 + a A a 5 (B 2 2 + Ci) 2 ](B 1 B 2 + dC 2 ) 

- ^a 2 a^B 1 B 2 + dC 2 + B X C 2 - B 2 Ci) 
25o 

x (B 1 B 2 + C ± C 2 - B X C 2 + B 2 C X ) [B\ + £? 2 + C\ + C 2 ) ( B 26) 

- ^[a ia3 (Bf - d){Bl - CD + a 3 a 5 (B 2 - C 2 )(5 2 4 - C 2 4 )] 

- l[a iaz {Bl + C\) + a 3 e* 5 (B 2 2 + Cf )]S 1 S 2 C 1 C 2 

o 

5 

- ^(«i« 4 + a 2 a 5 )(5ifi 2 + CxC,) 3 

- A {a 2 a 3 [(Bf + Clf + 3(S?S| + C?Cf) 

+ 2{B 1 C 2 + 5 2 Ci) 2 - 2fi 1 5 2 CiC 2 ] 
+a 3 «4[(5 2 2 + C 2 2 ) 2 + 3(BfB 2 2 + C 2 Cf) 

+ 2(B ± C 2 + 5 2 Ci) 2 - 2B 1 B 2 C 1 C 2 }} (B X B 2 + C X C 2 ). 



